Algebraic vortex liquid in spin- 1/2 triangular antiferromagnets: Scenario for CS2CUCI4 
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Motivated by inelastic neutron scattering data on CS2CUCI4, we explore spin-1/2 triangular lattice antiferro- 
magnets with both spatial and easy-plane exchange anisotropies, the latter due to an observed Dzyaloshinskii- 
Moriya interaction. Exploiting a duality mapping followed by a fermionization of the dual vortex degrees of 
freedom, we find a novel "critical" spin-liquid phase described in terms of Dirac fermions with an emergent 
global SU(4) symmetry minimally coupled to a non-compact U(l) gauge field. This "algebraic vortex liquid" 
supports gapless spin excitations and universal power-law correlations in the dynamical spin structure factor 
which are consistent with those observed in CS2CUCI4. We suggest future neutron scattering experiments that 
should help distinguish between the algebraic vortex liquid and other spin liquids and quantum critical points 
previously proposed in the context of CS2CUCI4. 
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The search for two-dimensional (2D) spin liquids has been 
one of the most tantalizing pursuits in condensed matter 
physics. Among the most promising systems for realizing 
such states is the spin-1/2 triangular antiferromagnet, as both 
the low spin and geometric frustration suppress magnetic or- 
dering. This was appreciated over three decades ago by An- 
derson, who postulated a quantum-disordered "resonating- 
valence-bond" (RVB) ground state in the spin-1/2 Heisen- 
berg triangular antiferromagnetil| . Anderson's RVB concept 
matured in the high-Tc era, with the triangular lattice of- 
ten center stage. Using slave bosons, Sachdev explored an 
Sp(Af) generalization of the Heisenberg antiferromagnet and 
obtained a quantum-disordered ground state, the Z2 spin liq- 
uid, which breaks no symmetries! 2]. More recently, Moess- 
ner and Sondhi realized a Z2 spin liquid in a quantum dimer 
model on the triangular lattice fs']. By exploiting a correspon- 
dence between the triangular antiferromagnet and hard-core 
bosons in a magnetic field, Kalmeyer and Laughlin introduced 
a "chiral" spin liquid|4] which violates time-reversal symme- 
try. Both the Z2 and chiral spin liquids admit gapped, frac- 
tionalized s = 1/2 excitations — spinons — which are bosonic 
in the former case and "semionic" in the latter. 

In spite of these theoretical advances, experimental spin- 
liquid candidates have only recently appeared. One promising 
system is the spin-1/2 triangular antiferromagnet CS2CUCI4, 
which has anisotropic exchange energies J = 4.3K and 
J' — 0.34 J (see Fig. 0. Although this material develops 
long-range spiral order below the Neel temperature T/v = 
0.62K|5], unusual features reminiscent of spin-liquid physics 
are manifested in its dynamics, probed via inelastic neutron 
scattering|5, 6]. Most notably, in addition to the sharp low- 
energy spin-wave peaks observed in the ordered phase, neu- 
tron scans at higher energies reveal "critical" power laws in 
the dynamical structure factor This enhanced scattering per- 
sists in a range of temperatures above T/v, where the magnons 
are absent, and is suggestive of spinon deconfinement charac- 
teristic of spin liquids. 

This remarkable behavior has attracted much theoretical in- 



terest in CS2CUCI4, and several scenarios for the origin of 
the anomalous scattering have been proposed. Spin-wave 
theorylJI 13, 0] and series expansion studies fiol fllll have 
yielded important quantitative connection with experiment. 
Quasi- ID effects have been explored by approaching the tri- 
angular lattice by coupling ID chains ildri3r . Sachdev's slave 
boson approach was generalized to the anisotropic triangular 
antiferromagnet by Chung et al. \ 14], and Isakov et al. \ 15] ex- 
plored the possibility that the CS2CUCI4 phenomenology may 
be controlled by a quantum critical point separating the Z2 
spin liquid and the spiral state. Using slave fermions, Zhou 
and Wen] 16] alternatively suggested the presence of a "criti- 
cal" algebraic spin liquid. 

Here, we pursue a new theoretical approach to the triangu- 
lar antiferromagnet. We consider an easy-plane XXZ spin-1/2 
system reformulated in terms of fermionized vortex degrees of 
freedom using Chern-Simons flux attachment. Remarkably, 
this approach leads naturally to a new "critical" spin liquid — 
the "algebraic vortex liquid" — which we explore and then ap- 
ply to CS2CUCI4. 

Algebraic Vortex Liquid. — Consider an easy-plane spin- 
1/2 antiferromagnet on the anisotropic triangular lattice shown 
in Fig.n We return below to the appropriateness of the easy- 
plane assumption for CS2CUCI4. We follow closely a dual ap- 
proach employing fermionized vortices, developed for integer- 
spin systems in Ref. 

IB 

Implementing the standard duality 
mapping, one obtains a theory for interacting bosonic vortices 
on the dual honeycomb lattice (see Fig.Q. A crucial feature 
is that the vortices are at half-filling due to the spin frustration. 
Moreover, because is a half-integer, the vortices "see" an 
average background of tt flux through each hexagon. In terms 
of a vortex number operator A'x and its conjugate phase ^x, 
the vortex Hamiltonian is 

Wdual = — ^ ^xx' C0s(6'x — ^x' — flxx' — ^xx') 
<xx'> 

+ 51 - l/2)14x'(iVx' - 1/2) + Ha ■ (1) 
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FIG. 1: Triangular lattice and the dual honeycomb on which vor- 
tices reside. Spins shown illustrate a vortex. The dominant spin 
coupling/vortex hopping occurs along the corresponding bold links. 
With our convention the vortices are at half-filling, e.g., in the spiral 
state the vortex number is one (zero) on up (down) triangles. 

Here T4;x' encodes a logarithmic vortex repulsion; a°^, is a 
static field satisfying (A x aP)^. = tt, where (A x aP)^ is 
a lattice curl around the hexagon encircling triangular lattice 
site r; and Tia describes the dynamics of the dual gauge field 
Oxx' residing on honeycomb links. The spin component 
appears here as a dual gauge flux, 5^ ~ (A x a)r/ (27r). The 
vortex hopping amplitudes txx' are anisotropic to reflect the 
lattice anisotropy of the exchanges. Thus, vortices hop more 
easily across weak spin links J' than across strong links J. 

Since the vortices are at half-filling, the dual theory appears 
as intractable as the original spin model. One can, however, 
make significant progress hyfennionizing the vortices. While 
similar approaches employing fermionized spins have been 
pursued! 18|, our treatment is appealing because vortex den- 
sity fluctuations are greatly suppressed by their logarithmic 
repulsion. Following Ref. 17, we perform flux attachment to 
convert the bosonic vortices into fermions coupled to a Chern- 
Simons gauge field. Taking a "flux-smeared" mean-field state, 
the Chern-Simons flux averages to 2tt per hexagon, which can 
be ignored on the lattice. Thus, the mean-field Hamiltonian 
describes half-filled fermionic vortices hopping on the hon- 
eycomb in a background of vr flux (due to a^^,). The cor- 
responding band structure reveals four gapless Dirac points. 
Focusing on low-energy excitations in the vicinity of the four 
nodes and including gauge fields, we obtain the Euclidean La- 
grangian density 

1 i 
+ TTji'^t.uxd^ax)^ + —ef^uxAfj^d^Ax + Af ■ (2) 
2e"^ 47r 

Equation (|2ji describes four flavors of two-component Dirac 
fermions ipa, a — 1,...,4, minimally coupled to a non- 
compact U(l) gauge field and a Chern-Simons field A^^. 
The gauge field mediates the logarithmic vortex repulsion, 
while the Chern-Simons terms implement the flux attachment. 
Here and below, the nodal velocity anisotropics are absorbed 
in the scaling of coordinates, while the Maxwell term is only 
schematic. Symmetries of the original XXZ spin model pre- 



clude all possible mass terms in Eq. 

The term €41 represents symmetry-allowed four-fermion 
terms arising from short-range vortex interactions. If suffi- 
ciently strong, such terms can drive spontaneous fermion mass 
generation, leading to various symmetry-broken states LI9I1. 
For example, giving the fermions a mass by introducing a 
staggered chemical potential on the honeycomb induces a 
vortex density wave coiTesponding to the spiral spin-ordered 
state. On the other hand, the state with a spontaneously gen- 
erated ip^tpa mass {i.e., all fermion masses having the same 
sign) corresponds to the Kalmeyer-Laughlin chiral spin liquid. 

In addition to broken-symmetry states, Eq. (|2ji allows us to 
access a new "critical" spin liquid. To this end, we rewrite the 
Lagrangian in terms of = + and integrate out the 
Chern-Simons field to obtain £ — £qed3 + 0{d'^cP'), with 

-Cqeds ='4^a{^~ i^)ipa + -^{^tiuxduaxY + A/- (3) 

Remarkably, up to higher-derivative terms which are hence- 
forth dropped, we arrive at non-compact quantum electrody- 
namics in 2h-1 dimensions (QED3), with = 4 flavors of 
two-component Dirac fermions. |22| Physically, vortex den- 
sity fluctuations are suppressed so strongly by interactions that 
exchange statistics play only a minor role. 

QED3 has been widely studied {e.g., see refs. in 1 17]), and 
in the large- limit realizes a nontrivial stable critical phase. 
For N < Nc, with some unknown Nc, it is believed that four- 
fermion terms become relevant and spontaneously generate 
fermion masses, destroying criticality (except at fine-tuned 
critical points). Here, we proceed with the assumption that 
Nc < 4, implying the presence of a stable critical phase for 
our dual fermionized-vortex theory. We now explore some of 
the properties of this "Algebraic Vortex Liquid" (AVL). 

Critical spin correlations in the AVL. — The AVL respects 
all symmetries of the microscopic spin system, exhibiting no 
magnetic or other types of order Since the Dirac fermions are 
gapless, power-law spin correlations are expected. Consider 
first the in-plane spin components. The spin raising opera- 
tor 5+ adds S^- = 1 and hence 27r dual gauge flux. Near 
this flux insertion the fermionic Hamiltonian has four zero- 
energy modes, one for each flavor. Physical (gauge-invariant) 
states are obtained by occupying two of these, so there are 
six distinct such "monopole insertions." Following the pro- 
cedure of Ref. we determine the momenta caiTied by the 
monopoles: two occur at the spiral ordering wave vectors ±Q 
and three occur at the midpoints Mi 2.3 of the Brillouin zone 
(BZ) edges (see Fig.|2}. Numerical diagonalization suggests 
that the sixth monopole, which carries zero momentum, does 
not have the same symmetry as .S+ HQfl . and thus by itself will 
not contribute to the in-plane spin correlations. 

Monopole insertions are known to have nontrivial power- 
law coiTelations in the large- A^ QED3 theory. The scaling 
dimension, A,n, of the monopoles can be estimated from the 
leading large-A^ result] 20|, 2 Am « 0.53A^. A naive extrap- 
olation to Af = 4 yields r/m = 2 Am - 1 « 1.12. Since aU 
monopoles have the same scaling dimension, S'^ is expected 
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FIG. 2: Momenta in the triangular BZ at which dominant power-law 
correlations in the dynamic spin structure factor occur in the 
AVL. We predict the same power law at all five momenta. For S^^, 
the dominant scattering is at ±Q and has a faster-decaying power 
law. Dashed lines show kx = ±7r characteristic of a ID system. 



to exhibit the same universal power-law correlations at all five 
momenta K.j shown in Fig.|2] Specifically, near each Kj , the 
dynamic spin structure factor is predicted to scale as 

5+-(k = K, + q, = Ak, (^fi"q,pll/, . (4) 

The amplitudes Ak are sensitive to short-distance physics 
and can differ significantly among the five wave vectors, par- 
ticularly in the anisotropic system. With J ^ J' the ampli- 
tude at wave vector M3 with fc^. = is expected to be much 
suppressed compared to the other four momenta, as the latter 
are near = tt where the dominant antiferromagnetic corre- 
lations occur along the nearly decoupled chains. 

The spin correlation behaves rather differently. At zero 
momentum, the correlation is that of the conserved dual gauge 
flux. However, a more prominent power law occurs at the spi- 
ral ordering wave vectors ±Q. These arise because an ex- 
pression for 5*^ in terms of continuum fields allows a term 
g'*^ Wa/3 V/3 +H.C. with a fermionic bilinear tpWip whose 
correlation is enhanced by gauge field fluctuations. Here, W 
can be obtained by considering a perturbation to the fermionic 
hopping Hamiltonian that adds static gauge flux through the 
hexagons modulated at wave vector Q. The iS^^(k, w) dy- 
namical spin structure factor near ±Q has a form similar to 
Eq. (|4}, except the corresponding anomalous dimension is that 
of an enhanced fermionic bilinear in QED3. Using the 
result|21], ?7cnh « 3- 128/(37r2A^), we estimate r/enh « 1-92. 
At other wave vectors the dynamical correlation (k, uj) ex- 
hibits subdominant power laws. 

CS2CUCI4 Hamiltonian. — To address the possible ex- 
perimental relevance of the AVL phase, we now consider 
the measured CS2CUCI4 Hamiltonian. Besides the Heisen- 
berg exchange with J = 4.3K, J' = 0.34 J, there is also a 
Dzyaloshinskii-Moriya (DM) interaction 



DM 



= — ^ D ■ Sr X (Sr+5i + Sr+52 



(5) 



The vectors 2 connect spins on neighboring chains as in 
Fig-E and D = Dz is oriented perpendicular to the triangular 
layers, with D — 0.053 J. Although there is a small interlayer 
coupling! 8|, we focus on the 2D system. 

Significantly, the DM term provides an easy-plane 
anisotropy, breaking the SU(2) spin symmetry of the Heisen- 
berg exchange down to U(l), and also violates inversion sym- 
metry r -r. Thus for T < Tn = 0.62K the DM term 
determines both the ordering plane, which coincides with the 
triangular layers, and the sign of the ordering wave vector, 
which along with D changes sign from one layer to the next. 
Despite the small value of D, the easy-plane anisotropy is am- 
plified since the DM interaction is not frustrated near the dom- 
inant antiferromagnetic wave vector along the chains (whereas 
the J' coupling is frustrated). To quantify this, we briefly con- 
sider a classical 2D spin system with CS2CUCI4 parameters. 
Without the DM term, a J — J' Heisenberg spin system re- 
mains disordered at all temperatures. With DM coupling the 
system has only U(l) spin symmetry and thus exhibits a low- 
temperature phase with quasi-long-range order (QLRO). We 
performed a classical Monte Carlo study and found this tran- 
sition at Tc « 0.27 JS^. A simple classical ground state anal- 
ysis indicates that about half of the "phase stiffness" originates 
from the DM coupling itself. Taking S*^ = 3/4 appropriate for 
spin- 1/2, we estimate that a single layer would obtain QLRO 
below Tc = 0.84K. Once each layer has QLRO, an arbitrarily 
small interlayer coupling would induce 3D long-range order, 
even though the interlayer coupling is frustrated in CS2CUCI4 
due to the alternating sign of _D|8]. This suggests that the ob- 
served spiral ordering at in CS2CUCI4 is primarily driven 
by the easy-plane character of spins in each layer. As the clas- 
sical treatment neglects quantum fluctuations, it is reasonable 
that the estimated Tc somewhat exceeds T^. These consider- 
ations suggest that vortices, required to drive the 2D ordering, 
acquire integrity as degrees of freedom in CS2CUCI4. 

Scenarios for AVL in CS2CUCI4. — To specifically address 
the applicability of the AVL to CS2CUCI4, we consider the ef- 
fect of the DM term on the AVL. In the fermionized-vortex 
Lagrangian the DM interaction appears as an inversion- 
breaking fermion mass term corresponding to a staggered vor- 
tex chemical potential, with bare mass m ~ D. This mass 
drives the system to the spiral state as indicated by the vertical 
flow in Fig. 13 with chirality dictated by the sign of D. Thus, 
with the DM interaction the observed spiral ground state of 
CS2CUCI4 emerges naturally out of the critical AVL. 

Since the DM term induces an easy-plane spin anisotropy 
and breaks inversion symmetry, applicability of the AVL to 
CS2CUCI4 requires a delicate balance. For the AVL to ap- 
ply on intermediate energy scales, the DM interaction must 
first produce sufficient easy-plane anisotropy for the descrip- 
tion in terms of vortices to be appropriate, before destabiliz- 
ing the AVL state toward the spiral order. This is scenario 1 
in the schematic flow diagram of Fig.|3l The alternative sce- 
nario 2 does not approach the easy-plane fixed point but is 
driven directly to the magnetic order In the latter case, in- 
termediate energy scales may be governed by an (unknown) 
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FIG. 3: Schematic renormalization group flows near the AVL fixed 
point in the presence of the DM interaction. The vertical axis is the 
inversion-breaking fermion mass m arising from the DM term. The 
bare CS2CUCI4 Hamiltonian is indicated with the star. For flow tra- 
jectory 1, the dynamics at intermediate energy scales is governed by 
the AVL theory. For trajectory 2, intermediate energy scales are con- 
trolled by an (unknown) SU(2)-invariant criticality. 



SU(2)-invariant criticality indicated with question marks in 
the figure. Below, we pursue the consequences of scenario 1 . 

Comparison with experiments. — The spin dynamics in 
CS2CUCI4 was measured in neutron scattering experiments by 
Coldea et alM\^- For T < Tn, CS2CUCI4 has well-defined 
spin waves, gapless near zero momentum and observed with a 
small gap near the ordering wave vectors ±Q, presumably due 
to weak violations of SU(2) spin symmetry by the DM term. 
Experiments also see fairly small spin-wave gaps at momenta 
Mi_2 of Fig. 12] (Series expansion studies of the Heisenberg 
model show that the gaps near these momenta and also M3 
deviate strongly from spin-wave theory 1 11].) Broad contin- 
uum scattering is observed above the spin-wave gaps near mo- 
menta ±Q and Mi 2, and persists even for T > T^. Notably, 
Ref. 6 reports power-law line shapes in scans near these mo- 
menta (scans J and G in Ref. 6), each with the same exponent, 
Vexp = 0.74. The AVL also admits gapless spin-1 excitations 
with power-law scaling at these momenta, with the leading es- 
timate of 77ni = 1.12 somewhat larger than experiment. These 
momentum space locations of enhanced scattering are deter- 
mined by physics on the scale of several lattice spacings, and 
provide important constraints on theory. The AVL theory cap- 
tures this aspect of the CS2CUCI4 phenomenology rather well. 

Since the above data are near — tt, which is the domi- 
nant wave vector in the quasi- ID limit, some caution is nec- 
essary. Specifically, with J « 3 J' strong contributions on in- 
termediate energy scales from antiferromagnetic correlations 
along the chains are not unexpected. A measurement midway 
between points Q and Mi would help determine the trans- 
verse /cy-dependence of the continuum scattering and clarify 
whether it is meaningful to speak of enhanced scattering near 
discrete momenta in the 2D BZ. The AVL also predicts en- 
hanced iS^ (k, oj) near M3, albeit with a smaller amplitude. 
Being less influenced by quasi- ID effects, further measure- 
ments near M3 would be useful. Polarized neutron exper- 
iments to search for the distinct easy-plane character of the 
AVL would also be very interesting. 

Competing 2D theoretical proposals with full SU(2) spin 
symmetry include an algebraic spin liquid ("UlC") proposed 



by Zhou and WenfT^ and the quantum critical point (QCP) 
scenario of Isakov et a/. jH^. Each makes distinct predictions 
for the momenta of the low-energy spin-1 excitations, which 
also differ from the AVL. Such characterizations can in princi- 
ple be used to discriminate among different theories, although 
the limited energy window for observing the continua and the 
material's strong anisotropy are unavoidable complications. 
Further experimental and theoretical studies should help clar- 
ify the true "spin-liquid" nature of this interesting material. 

Acknowledgments. — We would like to thank Leon Ba- 
lents, T. Senthil and Martin Veillette for sharing their insights, 
and especially Michael Hermele for an initial collaboration. 
This work was supported by the National Science Foundation 
(J. A.) through grants PHY-9907949 (O. I. M. and M. R A. F.) 
and DMR-0210790 (M. P A. F). 



[1] P. W. Anderson, Mater. Res. Bull. 8, 153 (1973). 
[2] S. Sachdev, Phys. Rev. B 45, 12377 (1992). 
[3] R. Moessner and S. L. Sondhi, Phys. Rev. Lett. 86, 1881 (2001). 
[4] V. Kalmeyer and R. B. Laughlin, Phys. Rev. Lett. 59, 2095 
(1987). 

[5] R. Coldea, D. A. Tennant, A. M. Tsvelik, and Z. Tylczynski, 

Phys. Rev. Lett. 86, 1335 (2001). 
[6] R. Coldea, D. A. Tennant, and Z. Tylczynski, Phys. Rev. B 68, 

134424 (2003). 
[7] A. E. Trumper, Phys. Rev. B 60, 2987 (1999). 
[8] M. Y. Veillette, J. T. Chalker, and R. Coldea, Phys. Rev. B 71, 

214426 (2005). 

[9] M. Y. Veillette, A. J. A. James, and F. H. L. Essler, cond- 
mat/0506667 (unpublished). 

[10] W. Zheng, R. H. McKenzie, and R. R. P Singh, Phys. Rev. B 
59, 14367 (1999); W. Zheng, R. R. P Singh, R. H. McKenzie, 
and R. Coldea, Phys. Rev. B 71, 134422 (2005). 

[11] W. Zheng, J. O. Fjaerestad, R. R. P Singh, R. H. McKenzie, 
and R. Coldea, cond-mat/0506400 (unpublished). 

[12] M. Bocquet, F. H. L. Essler, A. M. Tsvelik, and A. O. Gogolin, 
Phys. Rev. B 64, 094425 (2001). 

[13] O. A. Starykh and L. Balents (unpubhshed.). 

[14] C. H. Chung, J. B. Marston, and R. H. McKenzie, J. Phys.: 
Condens. Matter 13, 5159 (2001); C. H. Chung, K. Voelker, 
and Y. B. Kim, Phys. Rev. B 68, 094412 (2003). 

[15] S. V. Isakov, T. Senthil, and Y. B. Kim, cond-mat/0503241 (un- 
published). 

[16] Y. Zhou and X. G. Wen, cond-mat/02 10662 (unpublished). 
[17] J. Alicea, O. I. Motrunich, M. Hermele, and M. P. A. Fisher, 

Phys. Rev. B 72, 064407 (2005). 
[18] A. Lopez, A. G. Rojo, and E. Fradkin, Phys. Rev. B 49, 15139 

(1994); K. Yang, L. K. Warman, and S. M. Girvin, Phys. Rev. 

Lett. 70, 2641 (1993). 
[19] J. Alicea, O. I. Motrunich, and M. P. A. Fisher, in preparation. 
[20] V. Borokhov, A. Kapustin, and X. Wu, J. High Energy Phys. 11, 

049 (2002). 

[21] W. Rantner and X.-G. Wen, Phys. Rev. B 66, 144501 (2002). 

[22] Dropping the higher-derivative term is justified by power count- 
ing, but is potentially dangerous since it changes the formal 
symmetries of the continuum action. We assume that it is fine 
to do so when describing critical fluctuations. 



